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ABSTRACT 

We study the velocity distribution of Milky Way disk stars in a kiloparsec-sized region around 
the Sun, based on ~ 2 million M-type stars from DR7 of SDSS, which have newly re-calibrated 
absolute proper motions from combining SDSS positions with the USNO-B catalogue. We esti- 
mate photometric distances to all stars, accurate to ~ 20%, and combine them with the proper 
motions to derive tangential velocities for this kinematically unbiased sample of stars. Based on 
a statistical de-projection method we then derive the vertical profiles (to heights of Z = 800 pc 
above the disk plane) for the first and second moments of the three dimensional stellar velocity 
distribution. We find that (W) — —7 ± 1 km/s and {U) = — 9 ± 1 km/s, independent of height 
above the mid-plane, reflecting the Sun's motion with respect to the local standard of rest. In 
contrast, (V) changes distinctly from —20 ± 2 km/s in the mid-plane to {V) — —32 km/s at Z 
= 800 pc, reflecting an asymmetric drift of the stellar mean velocity that increases with height. 
All three components of the M-star velocity dispersion show a strong linear rise away from the 
mid-plane, most notably azz, which grows from 18 km/s {Z — 0) to 40 km/s (at Z = 800 pc). 
We determine the orientation of the velocity ellipsoid, and find a significant vertex deviation of 
20 to 25 degrees, which decreases only slightly to heights oi Z = 800 pc. Away from the mid- 
plane, our sample exhibits a remarkably large tilt of the velocity ellipsoid towards the Galactic 
plane, which reaches 20° at Z — 800 pc and which is not easily explained. Finally, we determine 
the ratio cr'^^/aj^n near the mid-plane, which in the epicyclic approximation implies an almost 
perfectly flat rotation curve at the Solar radius. 

Subject headings: Galaxy: kinematics and dynamics, (Galaxy:) solar neighbourhood 



1. Introduction 

The velocity distribution of stars in the Milky 
Way has been most comprehensively studied in 
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the solar vicinity. One of the most long standing 
systematic collection o f relevant da ta is the Cata- 
logue of Nearby Stars ( Glieseiri957 ). with a fourth 
edition {CNS4) now nearly completed (JahreiB, in 
preparation). The CNS4 contains space velocity 
data for a few thousand stars. A much larger 
data set, based on stars with proper motions mea- 
sured by the Hipparcos satellite (ESA 1997), was 
analyzed by Dehnen fc Binnev ( 19981 ). From the 
proper motions of ^ 10^ stars they could recon- 
struct the first and second moments of the ve- 
locity distribution using a statistical de-projection 
method which makes use of the many viewing di- 
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rections through the velocity distribution of the 
stars. The velocity dispersions of stars in the so- 
lar vic inity were also studied by Nord strom et all 
(j2004l ) from the Geneva- Copenhagen Survey of F 
and G stars. These authors complemented the 
parallaxes and proper motions of the stars in the 
Hipparcos catalogue by their own radial veloc- 
ity measurements. With the advent of the deep, 
wide-field SDSS and RAVE surveys it has be- 
come possible to study the kinematics of stars in 
the Galactic disk a t much larger di s tances from 
the Sun. Recently iBochanski et al. ( 2007 ) have 
drawn from the Fifth Data Release of the Sloan 



Digit ized Sky Survey (jAdelman-McCarthv et al 



20071 ) a sample of about 7000 K and M stars 
for which three-dimensional space velocities are 
available (the SLoMaSS sample). These stars 
were selected along a line-of-sight which reaches 
a distance of 1 kp c ab ove the Galactic midp lane. 
Veltz et al.l (|2008l ) and lKlement et al.l (|2008f) have 



anal yzed the First Data Re l ease of the RAVE S ur- 
vey Isteinmetz et al.ll2006l ). IVeltz et al.l (|2008l ) in- 
ferred the kinematics of stars in the solar cylin- 
der by modeling the distribution of stars in phase 
space, which was t hen constrained by st ar counts 
of the RAVE data,. iKlement et"all (|2008l) searched 
for fin e structure in local phase space. Most re- 
cently, [SieberFetin- (2Qq3) have analyzed the ve- 
locity distribution of red clump giants at heights 
between 500 and 1500 pc above the Galactic mid- 
plane, using data f rom the Second Dat a Release of 
the RAVE Survey (jZwitter et al.ll2008l ). 

In this paper we analyze stellar data from the 
Seventh Data Release of the Sloan Digital Sky 
Survey (SDSS) (jAbazajian et al.ll2008l ). which in- 
cludes imaging taken as part of the Sloan Ex- 
tension for Galactic Under standing and Explo - 



( SEGUE) program (lYannv et all l2008h . 



ration 

one of three sub-surveys that were part of first 
extension of the SDSS, known as SDSS-II. The 
SEGUE program was designed, in part, to ob- 
tain ugriz imaging of some 3500 square degrees 
of sky outside of the original SDSS-I footprint 
( Fukugita eti r 'l996l Gunn et al . 1998 , 2006 



(iFukugita et a l.' 'l996l Gunn et al . 1998, 2006. 
York et al.l ijoO O, Hogg et al.l I2OOIL ISmith et al 
200l lStoughton el al.ll2qol A b azajianet al . 2003 . 
2004. 2005. IPier et al.l l2003l llvezic et al.l 12004 



Adelman-McCarth y et al. 2006, 2007, 2008, 
Tucker et al.1l2006l . The regions of sky targeted are 
primarily at lower Galactic latitudes <^ 35°), 



in order to better sample the disk/halo inter- 
face of the Milky Way. SEGUE also obtained 
R ~ 2000 spectroscopy, over the wavelength range 
3800 - 9200 A, for some 250,000 stars in 200 se- 
lected areas over t he sky available fr om Apache 
Point, New Mexico (jYannv et al.ll2008l ). The spec- 
troscopic candidates are selected on the basis of 
ugriz photometry to populate target categories 
chosen to explore the nature of the Galactic stel- 
lar populations at distances from 0.5 kpc to over 
100 kpc from the Sun. The SEGUE Stellar Pa- 
rameter Pipeline (SSPP; see Lee et al. 2008a, b. 



AUende Prieto et al.l l2008f ) obtained estimates of 



radial velocities, and in particular, the stellar at- 
mospheric parameters, for the targeted stars over 
wide ranges of effective temperature (Teg), surface 
gravity (log g), and metallicity ([Fe/H]). 

Of central importance to the present paper, the 
7 catalog contains astrometric and photomet- 
ric (ugriz) data for well over 100 million stars. 
Many of these have been iden tified in the USNO - 
B catalog of proper motions (|Monet et al. I l2003h . 
These measurements have been improved through 
a re-calibration effort based on SDSS data, and 
provi de us with an exte nsive set of kinematical 
data ( Munn et al.l |2004[ ). After correction of a 



small systematic error in the data reduction pro- 
cedure (Munn 2007, internal SDSS memorandum; 
note that the DR-7 CA^ proper motions have been 
corrected for this error), the absolute accuracy of 
the proper motions has reached now a few milliarc- 
seconds per year, corresponding to, e.g., ~5 km/s 
at a distance of 530 pc for — 2 mas/year. They 
are thus superbly suited to study the kinematics 
of the Galaxy, especially for stars that populate a 
local volume. 

We have extracted information for about 2 mil- 
lion M-type stars in the solar cylinder from the 
SDSS DiZ- 7 database; these stars represent a fair 
sample of an old, relaxed stellar population mostly 
drawn from the Galactic thin disk. Using a pho- 
tometric distance estimator we convert the proper 
motions into tangential velocities. For individ- 
ual stars, tangential velocities alone do not al- 
low the construction of space velocity components 
beca use the radial velo c ities a re missing. How- 
ever, Dehnen fc BinnevI ( 19981 ) have shown that, 
by making use of the many viewing angles through 
the velocity distribution of stars over a large frac- 
tion of the sky, the characteristic moments of the 
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velocity distribution can be determined in a statis- 
tical manner. Following this method, we are able 
to determine empirically, for the first time with 
high precision, the first and second moments of 
the velocity distribution of local stars as a func- 
tion of height above the Galactic midplane. 

Our paper is organized as follows. In §2 we de- 
scribe the sample selection and our methods for 
obtaining distances and deriving tangential veloc- 
ities. In §3 we apply th e statistical approach of 
Dehnen &: BinnevI ( 19981 ) to reconstruct the mo- 



( Abazaiian et al. 2005f) and determined spectro- 
sco pic parall axes of these stars. If we use fits to 



ments of the three-dimensional velocity distribu- 
tion of stars from their tangential velocities. Our 
conclusions are summarized in §4. 

2. Sample Selection, Distance Estimates, 
and Tangential Velocities 

In order to obtain a sample of stars that are 
representative of the old, relaxed population dom- 
inated by the thin disk of the Milky Way, we have 
extracted from the SDSS DR-7 database all ob- 
jects classified as stars in the de-reddened color 
range 0.49 < r-i < 1.6, 1.25 < g-r < 1.50, 
and with apparent in agnitudes 15 < g < 20.5 
fcf. Ilvezic et"aIll2Q05l) . These stars have spectral 
types MO and later. The upper cutoff a.t g — 20.5 
mag, which corresponds to r = 19.5 mag, is nec- 
essary because the proper motions are seriously 
contaminated by misidentifications at fainter mag- 
nitudes (Lepine 2008, internal SDSS memoran- 
dum). We have selected only stars with non-zero 
proper motion errors, apparent magnitude errors 
< 0.05 mag in the g,r,i bands, and which were 
not flagged for any known astrometric or photo- 
metric peculiarities. This resulted in a sample of 
about 2.6 million stars (mean age expected to be 
~5 Gyr). 

We have determined a photometric distance es- 
timate for each sta r using the absolut e magnitude- 
color relations of Juric et al. ( 2008h . We have 
compared both the 'faint' and 'bright' nor- 
malizations with color -mag nitude relat i ons b y 
Davenp ort et all (l2006l) and IWest et all (|2005l ). 



Davenport et al 



the IWest et a l. (2005.) or iDavenport et all (|2006h 



20061 ) have observed M stars 



data as distance estimators in the de-projection 
formalism (described below) to obtain the first 
and second moments of the velocity distribution 
of the stars, we find that the results are practi- 
call y indistingu i shable from the results based on 
the lJuric et all (|2008l ) bright normalization. We 
conclude that this normalization is thus well suited 
to serve as a photometric distance estimator for 
red main sequence stars in the Galactic thin disk, 
such as M dwarfs, 

Mr = 3.2 + 13.30 (r-i)- 11.50 (r-i)2 

+5A0{r -if -0.70 (r-i)'^ . (1) 

As explained below, we count M dwarfs up to dis- 
tances reaching 800 pc above the Galactic mid- 
plane. According t o the model of the G alactic thin 
and thick disks bv lJuric et al. (2008), we expect 
about 13% of the M dwarfs close to the midplane 
to be thick-disk stars; in the distance range of 600 
to 800 pc above the midplane the fraction of thick 
disk stars, which are more metal-poor than thin 
disk stars, is expected to be about 36%. There- 
fore, we have compared t h e fain t and bright nor- 
malizations of Juric et al. I dlooi) with the fiducial 
line delineating the main seq uence of the globu- 
lar cluster M71 as observed bv lClem et all (|2008l ). 
transformed into the native uqriz system using the 
relations of Tucker et al. (2006.), and assuming a 
distance modulus of (m — M)n = 13.02. M7 1 has a 
metalhcity of [Fe/H] = -0.73 (iHarris Ill996l ) which 
is typical for the majority of thick-disk stars. This 
comparison shows that for metal-poor stars in the 
color range r — i > 0.5, photometric distances can 
be reliably determined using the faint normaliza- 
tion, 

Mr = 4.0 + 11.86 (r-i)- 10.74 (r-i)2 

+5.99 (r-i)3- 1.20 (r-i)-*. (2) 

In the following we provide results using both dis- 
tance scales in order to demonstrate how (little) 



drawn from the Catalogue of Nearby Stars (jJahreifi &: 
19971 ). for which trigonometric par allaxes are 



Wielenf" ''^^^^^^ systematically affected by the pres- 
5nce of thick-disk stars. The difference between 



availa ble, in the ugriz filter system. IWest et al 



(|2005l ) have selected M dwarfs from the Third 
Data Release of the Sloan Digital Sky Survey 



the two calibrations is illustrated in Fig.[l]as func- 
tion of de-reddened (r — i)o color. 

The right ascension and declination, taken to- 
gether with the distance estimate, were used to 
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calculate the position of each star in a Cartesian 
reference system centered on the position of the 
Sun, with the X-axis pointing toward the Galactic 
center, the F-axis in the direction of Galactic ro- 
tation, and the Z-eods towards the Galactic north 
pole, respectively. We then defined eight pill-box 
shaped counting volumes, oriented parallel to the 
Galactic midplane and stacked onto each other. 
Each volume has a vertical height of 100 pc and a 
radial size \/X^ + Y"^ < 1 kpc. The midpoints of 
the counting volumes are centered on the Galacto- 
centric radius of the Sun and Z — 50, 150, ... 750 
pc, respectively. In Fig. [T] we have overlaid fre- 
quency distributions of the colors of stars in the 
distance bin closest to the Galactic midplane and 
in the most distant bin, respectively, on the abso- 
lute magnitude differences between the faint and 
bright calibrations. These indicate which color 
ranges of the calibration of the absolute magni- 
tudes are relevant in the various distance bins. 

Fig. [2] shows the frequency distributions of the 
proper motions in right ascension and declination 
for the stars in each counting volume. Averaged 
over the entire sample the errors of the proper mo- 
tions are e^^ and e^^ = 3.36 mas/year. 

Fig. [3] illustrates our data as a Hess diagram of 
the reduced proper motions for our stars 



Hr = r + 5 log fi + 5 



(3) 



versus their (r — i)o colors. In eq. [3] /i denotes 
the total proper motion of each star in units of 
mas/year. Even though the stars have been se- 
lected into the counting volumes according to their 
estimated distances, their reduced proper motions 



are t o first order distance independent (jLuvten 



19391) . and are thus also independent from the cho- 
sen distance scale. Since M giants are intrinsically 
brighter (by at least 9 magnitudes) than M dwarfs, 
they would show up in Fig. [3] in the lower flank of 
the distribution. As can be seen from Fig. [31 the 
contamination of our sample by su ch halo M giants 
is, as expected (|Covev et al.ll2008l ). negligible. 

Fig. [4] shows the distribution of the selected 
stars as function of height above the Galactic mid- 
plane versus their (r — i)o colors. As can be seen in 
this figure, our sample is volume limited for stars 
bluer than (r — i)o < 0.95 (with the exception of 
perhaps a few bright, very nearby stars). In the 
upper right portion of Fig. |31 red M stars are miss- 
ing whose apparent magnitudes are fainter than 



the magnitude limit of the survey. 

Returning to Fig. [3l wc note a mild Malmquist 
bias in our sample. The high altitude counting vol- 
umes are populated by comparatively more blue 
stars than red stars, and, given the apparent mag- 
nitude limits of the survey, the blue stars are on 
the average intrinsically brighter. 

For each star we have calculated its tangential 
velocity from the proper motions and the photo- 
metric distances b a sed o n the bright normaliza- 
tion bv lJuric et al. ( 2008 ). Assuming zero line-of- 
sight velocities, we then computed three- dimcn- 
sional space velocities. F ollowing the notation of 



Dehnen fc BinnevI (|1998f ). these are referred to as 
the proper motion velocities pu, Pv, and pw- The 
velocity components are oriented in the same way 
as the spatial coordinates. Finally, we rejected 
stars which had at least one single velocity com- 
ponent > 150 km/s. After application of this cut, 
our final sample comprises 1 898 598 stars. Their 
spatial distribution is shown in Fig. [5l The errors 
of the space velocity components were determined 
by application of a standard Gaussian error prot 



agation formalism. Following lJuric et al.l ()200 



we have assumed individual relative distance er- 
rors of 20%, which corresponds to an error in the 
absolute magnitudes of 0.43 mag. These errors re- 
sult in typical errors for the individual velocities of 
eptj' <^pv> and ep^, = 2.5{+2, -1) km/s, 2.5(-|-3, -1) 
km/s, and 1.5(-|-3,-l) km/s respectively, but with 
tails out to ^ 10 km/s. However, the individual 
errors are greatly diminished when one considers 
the large number of stars in each of our counting 
volumes. For instance, if we were to change the 
individual distance errors from 20 percent to 30 
percent, the resulting first and second moments 
of the velocity distribution would be precisely the 
same as given in Table 1. 

The last step in the preparation of the data was 
the correction of each individual proper motion 
velocity for Galactic mean rotation 

Pu = pu — B)d cos (6) sin (l) 

pv = pv~ {A + B)d cos (6) cos {I), (4) 

where d denotes the distance of a star and I 
and b are Galactic longitude and latitude, respec- 
tively. In this way we have constructed proper 
motion velocities that reflect the true peculiar 
velocities of the stars. For the Oort constants 
A and B we adopt the values determined by 
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Feast fc Whitelockl (|l997tl . A 
and B = -12.37 km/s/kpc. 



14.82 km/s/kpc 



3. Determination of the Velocity Distribu- 
tion Moments 

3.1. Projection Formalism 

The proper motion velocities, which we have 
introduced in the previous section, are formally 
defined as 



((C/,F,M^),e,) e„ (5) 



where eV is the unit vector pointing to the star 
tinder consideration, 



PU ^ 




( V 


Pv 




V 


Pw J 




\ w 



COS {h) cos (/) 
fV = I COS (6) sin (Z) 
sin {h) 



(6) 
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Equation [5] can be written as a set of algebraic 
equations 



(7) 



The elements of the symmetric matrix A are given 
in Appendix A. The key assumption of the de- 
projection of the proper motion velocities is that 
the lines of sight towards the stars are statisti- 
cally uncorrelated with t he velocities of the stars 
( Dehnen fc Binnev Il998l ). We thus take the sta- 
tistical average of eq. [7] over the samples in each 
counting volume, as defined above (-^ ^^^=1) 




{A) 




(8) 



In contrast to eq. [7l which cannot be inverted, 
eq. [8] can be inverted as 





(9) 



which gives the three first moments of the velocity 
distribution. Since we presume to have selected a 
sample from a relaxed population of stars in the 
Galactic disk, ([/) and {W) should reflect simply 



the solar motion relative to the standard of the 
rest, Uq = -{U) and VFq = -{W). The quantity 
{V) is, due to the asymmetry of the V velocity 
distribution, more negative than —Vq. 

Next we consider the second moments of the ve- 
locity distribution. The components of the proper 
motion velocity of each star can be combined as 



PiPj = ^AikVkAjmVm (10) 
/c,m 

^ ^^{AikAjni ^~ AjkAira)V}iVra^ 
k.m 

where the indices stand for i,j,k,m = U,V,W 
and {vu, vy, vw) = {U, V, W). Again we take the 
ensemble average 



1 



k.m 



(11) 



which is manifestly symmetric in the indices 
Formally, we set up a set of algebraic equations 



PuPu 
PuPv 
PuPw 
PvPv 
PvPw 
PwPw 



(B) 




(12) 



Since the lines of sight to the stars are statistically 
uncorrelated with their velocities, the matrix (B) 
is also not singular, and allows for each Z-slicc the 
inversion 




PuPu 

PUPV 

PuPw 
PvPv 

PVPW 

pwpw 



(13) 



The elements of matrix B are also given in Ap- 
pendix A. 

The velocity dispersions are given by crj^j^ — 
{UU)-{U)^ = {VV)~{V)\ 4z = (W^W^)- 
(VF)^ and we have the mixed moments = 
-(([/y)-(C/)(F)),4^ = -((t7l^)-(C/)(T^)),a2^ 
lyW) - (V^)(W^). In the definition of the latter, 
we have switched from U to R — ~U which is, for 
instance, more convenient to use with the Jeans 
equations. 
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It seems intuitive that this inversion should 
work for an all-sky sample, but less so in the case 
of the SDSS coverage. In order to test that the 
sky coverage of the selected sample of stars is suf- 
ficient to reliably invert eqns. [5] and [T^l we have 
carried out a Monte-Carlo simulation. We have as- 
signed to each star in the counting volumes a space 
velocity drawn randomly from the Schwarzschild 
distribution 



/ oc exp - 



u-{u) 



v-{v)V fW-{W) 



(14) 



(^4,4, 



with ([/) ^ -10 km/s, (y) = -26 km/s, (VF> = -7 
km/s, a^,ji — 45 km/s, cr^^ = 32 km/s, and azz = 
24 km/s. We have converted the space velocities to 
proper motion velocities and found that we could 
reconstruct, from just the samples contained in the 
SDSS footprint, both the assumed velocity offsets 
and velocity dispersions to accuracies of 0.2 km/s, 
respectively. 

3.2. Results 

The de-projection procedure provides us with 
an observational estimate of the stellar velocity 
dispersion ellipsoid as a function of height above 
the Galactic midplane. Results are given in Table 
1 and illustrated in Fig. [B] The error calculus for 
all entries in the table is described in Appendix B. 
As can be seen in Fig. [SI there is a mild scatter 
of the data points which we ascribe to individual 
distance errors. 

In order to demonstrate the systematic effect of 
the adopted distance scale, we also show in Fig. [S] 
kinematical results based on t he faint nor maliza- 
tion of the main sequence by lJuric et al. I Hooa). 
As discussed in §2, this distance scale is appro- 
priate for metal-poor thick-disk stars. Since our 
sample is presumed to be significantly contami- 
nated by thick-disk stars only in the most distant 
counting volumes, the alternative distance scale 
should be actually used only in these counting vol- 
umes. As can be seen from the panels on the right 
hand side of Fig. [6l switching the distance scale 
leads to differences of the the mean motions and 
velocity dispersions of 1 to 2 km/s, with the ex- 
ception of CTflH, which deviates by up to 3 km/s. 



These changes are unexpectedly small, but can be 
understood from closer inspection of Fig. [1] The 
lowest distance bin is populated by comparatively 
red stars with typical colors 1.1 <^ {r — i)o <J 
<1.5. As can be seen from Fig. [1] the difference 
Mr(faint) — Afr(bright) ranges in this color range 
from 0.2 to -0.2, mag so that in the statistical aver- 
age no change of the kinematic results is expected. 
In the most distant bin, 0.5 < {r — i)o < 1, the dif- 
ference Mr(faint) — Mr(bright) is about 0.2 mag, 
which corresponds to a distance scale contraction 
of 9%. However, as a result of the reduced distance 
estimates, stars will wander into the Z-slice from 
above, and some stars will drop through the bot- 
tom of the counting volume. As can be seen from 
Fig- HI the incoming stars (from above) will be sys- 
tematically bluer than the stars which they replace 
(from below). The former will also be systemat- 
ically brighter, which compensates to a consider- 
able degree the switch from the bright to the faint 
calibration of the absolute magnitudes. Thus, the 
systematic errors of the kinematical data intro- 
duced by systematic errors in the distance scale 
are only of order ~ 1 to 2 km/s. 

As an independent consistency check on our 
results, we have retrieved from the SpecPhoto 
database of SDSS DR- 7 the radial velocities of all 
stars at galactic latitudes b > 70° in the color win- 
dow 1-0 < g — r < 1.5 and 0.3 < r — z < 1.6. Since 
these stars lie close to the Galactic pole, their ra- 
dial velocities reflect their W velocity components. 
The resulting mean velocities and velocity disper- 
sions are summarized as a function of height above 
the Galactic midplane in Table 2. As can be seen 
from a comparison with Table 1 , the vertical veloc- 
ity dispersions determined in this way are in excel- 
lent agreement with those estimated from proper 
motions alone. 

The second row in Table 1 lists, for compari- 
son, kinematical data for late-type star^ repro- 
duced from the CNS4 (JahreiB, in preparation), 
which contains data for stars within a distance 
of 25 pc from the Sun. With the exception of 
the mean radial velocity ([/), the other first and 
second moments of the velocity distribution de- 
termined in the Z = 0-100 pc sample of SDSS 
stars are remarkably consistent with the relatively 
tiny sample of nearby stars. The third row in Ta- 



1 Groups 2-5 as defined in ljahreifi fc WielenI l|l997f ) 
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Table 1 

First and second moments of the velocity distribution as function of height above the 

Galactic midplane 



z 

pc 


N 


km/s 


{V) 

km/s 


{W) 
km/s 


O'RR 

km/s 


km/s 


km,/ s 


{UV) 
(km/s)^ 


{UW) 
(km/ s)"^ 


{VW) 
(km/ s)"^ 


U-IUU 


A O ATC 

48 075 


-8.52 


OA A A 

-20.04 


^ 1 A 

-7. 10 


32.4 


23. U 


18.1 


OOA 'T 


32.1 


127.4 






n 99 


n 1 zi 

U. ±^4 


u. ±0 


nil 

U.l± 


n HQ 
u.uy 


n IS 

U.lo 




0. / 






oUu 


'12.2 


-22.3 


-0.9 


38.2 


25.0 


-in 1 
19.4 


151.1 


-2b.7 


-1.8 






9 1 


1 ^ 

1 . D 


1 ^ 


1 7 


1 1 
1.4 


1 . 












-1044 


-24.59 


-8.17 


00.8 


2o. 7 


18.3 


1 i /? n 

410.0 


-Ho.l 


202.7 






1 fiQ 


1 '71 
1.11 


1.10 


1.10 


1. IZ 


9 KR 
/O. Do 








1 r\r\ OAA 

100-21)0 


OOA AO A 

280 929 


O 1 A 

-8.10 


-21.62 


-7.14 


35.3 


24.9 


OA A 

20.9 


OAO A 

298.9 


7.1 


t AO 

143.3 






n no 


u.uo 


n ns 

U.UO 




n fizL 


n ns 

U.UO 


1 7 


1 7 


1 9 
l.Z 


200-300 


387 749 


-8.05 


-22.37 


-7.49 


38.5 


27.2 


24.9 


292.9 


-42.2 


149.6 






0.08 


0.05 


0.08 


0.04 


0.03 


0.08 


1.6 


1.7 


1.4 


300-400 


360829 


-9.06 


-23.51 


-7.29 


41.5 


29.8 


28.0 


325.5 


-75.9 


137.0 






0.08 


0.06 


0.09 


0.04 


0.04 


0.10 


1.9 


2.1 


1.9 


400-500 


302 601 


-9.62 


-24.76 


-7.35 


43.1 


31.7 


30.6 


332.2 


-125.4 


135.6 






0.09 


0.07 


0.11 


0.05 


0.05 


0.13 


2.3 


2.7 


2.5 


500-600 


237746 


-10.28 


-26.22 


-7.12 


45.0 


33.6 


33.3 


342.8 


-181.1 


122.3 






0.11 


0.09 


0.14 


0.06 


0.06 


0.17 


3.0 


3.6 


3.6 


600-700 


168 595 


-10.19 


-28.61 


-7.12 


47.0 


35.4 


36.6 


358.6 


-232.5 


132.3 






0.13 


0.11 


0.19 


0.08 


0.07 


0.23 


3.9 


4.7 


5.0 


700-800 


112074 


-10.42 


-32.18 


-7.18 


49.1 


37.6 


40.0 


419.2 


-263.0 


166.8 






0.17 


0.14 


0.27 


0.11 


0.09 


0.36 


5.4 


6.9 


7.7 



Col. (1): distance interval, col. (2): number of stars, cols. (3) to (5): mean velocity components and their 
errors (below), cols. (6) to (8): velocity dispersions and their errors (below), cols. (9) to (11): mixed 
moments of the velocity ellipsoid and their errors (below). The errors represent only internal errors. 

Systematic errors can be as large as 2-3 km/s. 



Table 2 

First and second moments of the distribution of spectroscopically measured line-of-sight 

velocities^ of m dwarfs at latitudes b > 70° 



z 


N 


{Vlos) 








pc 




km/s 


km/s 


km/s 


km/s 


0-200 


272 


-9.5 


0.6 


24.0 


1.0 


200-400 


770 


-12.0 


0.4 


27.6 


0.7 


400 600 


247 


-8.5 


0.5 


29.7 


1.3 


600-800 


64 


-16.0 


2.0 


39.6 


3.5 



tSpecPhoto data from SDSS DR-7 
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ble 1 gives the results of iDehnen fc Binnev (|l998l) 
in their reddest bin, which corresponds to K-type 
stars. Again, there is good agreement between 
their results and the results i n our low l a titude 
sample. We note as well that IVeltz et al. I (I2OO8I) 
find, from their analysis of the RAVE survey 
(fsteinmetz et al.ll2006h (U) = -8.5 km/s, which is 
very close to our SDSS result. Moreover, we can 
compare our result s with the obseryations of the 
SLoMaSS stars by iBochanski et al.l (|2007t ). For 
this purpose we have combined all 1 380 183 stars 
with < Z < 500 pc and derived with the sta- 
tistical de-projection method velocity dispersions 
i(TRR,cr^^,crzz) = (40.0, 28.5, 25.4) km/s. These 
ar e in excellent agreement with the measurements 
bv IBochanski et aL IIIOO^) (39.0, 30.0, 24.8) km/s. 

The distribution of the [/-velocity components 
of the stars is expected to be symmetrical with 
respect to the local standard of rest (LSR). Thus 
the mean velocity (U) reflects the solar motion 
relative to the LSR, (U) = —Uq. But as can be 
seen from Table 1, we find a slight, but statistically 
significant, decrease of (U) with increasing height 
above the midplane. This might indicate some 
radially inward streaming of the local stars relative 
to the rest of the disk, perhaps induced by a spiral 
arm. 

The mean rotation velocity (V) does not sim- 
ply reflect the solar motion — V0, but is much 
more negative due to the asymmetrical shape of 
the distribution of the IZ-velocity components of 
the stars. The solar motion in V alone would be 
reflected by a subsample of yo ung stars with very 
small velocity dispersions (cf. IPehnen fc Binnev 
19981 ). which we cannot isolate in our data, be- 
cause we do not know the individual space ve- 
locities of the stars. The fourth column of Table 
1 shows that the mean rotation velocity of stars 
slows down systematically with height above the 
midplane, as A{V)/Aa%j^ = -0.0089 (km/s)-\ 
and A{V)/AZ ^ - 0.017 km/s/kpc. The lat- 
ter is expected due to the asymmetric drift effect 
( Binnev fc Tremaind 1987h . A quantitative inter- 
pretation of this effect is not possible on the ba- 
sis of the present data. For Jeans modeling of 
the Milky Way disk, the radial gradient of the 
c_Rfl velocity dispersion must be known. Unfortu- 
nately, this cannot be measured using stars from 
the SDSS DR-7 sample, but we hope to address 
this problem in the near future. 



The quantity (W) = —Wq stays constant with 
height above the midplane. 

The velocity dispersions (trr, cr^^, and azz 
rise systematically with height above the mid- 
plane, Z. This rise is expected due to two effects. 
Even the relaxed population of M-type stars rep- 
resents an a ge mix of young and old stars. Since 
the work of Spitzer fc Sch warzschildl (1951), it is 
well known that older stars exhibit, as a result of 
'disk heating', significantly l arger velocity disp er- 
sions than younger stars (see iFuchs et al.ll2001 for 
a review). Stars with larger vertical azz veloc- 
ity dispersions can ascend higher into the Galac- 
tic gravitational potential, so that one expects 
to find stars with higher velocity dispersions at 
greater heights above the midplane than closer to 
the midplane (|Fuchs fc WielenI ll987L IWest et aL 
12008 ) . Furthermore, our sample transitions with 
increasing height above the midplane from the old 
thin-disk to the thick- disk population. Thick- 
disk stars have velocity dispersions on the order 
( 45-65, 40-50, 30-40) km/s 



of ((TRR,aM,crzz) = ( 
(IChiba fc Beers ' '2000", 



Alco be fc Cubarsi I I2005L 



Holmbcrg et al. 2 007 ), which are only slightly 
larger than what we actually measure. We have 
estimated, from the 50% contour level of the ve- 
locity distribution of red clump giants at alti- 
tudes between 500 p c and 1500 p c abov e the 
midplane derived by ISiebert et al.l ( 2008f l with 
RAVE DR-2 data, velocity dispersions of about 
{<^RR,<^4>4>,<^zz) ~ (79, 56, 35) km/s. The vertical 
velocity dispersion is consistent with our mea- 
surement in the highest altitude bin, whereas the 
planar velocity dispersions are significantly higher 
than our measurements, and, particularly in the U 
velocity component, even higher than what is com- 
monly ascribed to thick-disk stars. We note, how- 
ever, that such high ajiji and a^^ dispersions are 
co mparable with th e veloc ity dispersions ascribed 
bv IBochanski et aL ( 2007t ) to the high-dispersion 
wings of the velocity distribution o f the stars in 
their sample. Thus, the sample of ISiebert et al. 
(|2008l ). with a 'ceihng' at Z = 1500 pc, seems 



to be co ntaminated by hot metal-weak thick disk 
(|Chiba fc Beers 2000) and halo stars; the tilt of 
the velocity ellipsoid which they measure might 
not refiect that of disk stars. 

The rotational lag of the stars is tied by the 
asymmetric drift effect to the radial velocity dis- 
persion uuR, so that —(y) increases with Z, pre- 
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cisely as observed. As explained above, a quanti- 
tative modeling of this effect is not possible on the 
basis of the present data alone. 

Owing to the richness of our data, we were able 
to split them up into subsets based on color. As 
a test, we have selected only those stars with col- 
ors in the range 0.49 < (r — i)o < 0.95, where our 
sample is volume complete (cf. Fig. 2]). If we use 
only these stars in the de-projection of the proper 
motion velocities, wc find the first and second mo- 
ments of the velocity distribution illustrated in 
Fig. [71 A comparison with Fig. [6] shows that the 
kinematical results inferred from the 'blue' subset 
are entirely consistent with that of the total data 
set. Only in the (V) velocity component do we 
find, close to the midplane, an offset of 4 km/s. 

Next we discuss the off-diagonal elements of the 
velocity dispersion tensor. The mixed moment 
aR^ is illustrated in Fig. [8] as a function of height 
above the Galactic midplane. The inset of Fig. [5] 
shows the implied vertex deviation 



1 

- — arctan ■ 



(15) 



'RR 



As can be seen from Fig. |51 the vertex devia- 
tion is non-negligible, even if it falls slightly with 
height above the midplane. Close to the midplane 
our measurement i s cons istent with the result of 
Dehnen fc BinnevI (jl998[ ). Since our sample of 



stars is supposed to be drawn from a dynamically 
relaxed population of the Milky Way, the vertex 
deviation is almost certainly due to dynamical ef- 
fects induced by non-axisymmetric structures in 
the Galactic disk. 

The moment can be used to determine the 
apparent tilting angle of the velocity ellipsoid in 
the meridional plane spanned by U and W, 



1 



■ arctan ■ 



2 4z 



(16) 



'RR 



' ZZ 



The mixed moment itself is shown as — y crjj^ in 
Fig. [9] as a function of height above the Galac- 
tic midplane. The resulting tilting angle a is de- 
picted in the inset of Fig. [51 Close to the midplane 
we find a ~ implying that the velocity ellipsoid 
points towards the Galactic Center. Above the 
midplane the velocity ellipsoid is tilted towards 
the Galactic center, however, by a much larger an- 
gle than previously thought ( Binnev fc Tremaind 



1987LlGilmore1ll989llKent fc de Zeeuw Ill99lfl . At 

a height oi Z = 750 pc the velocity ellipsoid is 
tilted by a = —20°, whereas one would expect a = 
—5.4° if the velocity ellipsoid was pointing straight 
towards the Galactic center. In order to assess 
the plausibility of this unexpected result we ad- 
dress each of the three input parameters of eq. [TBj 
in turn. We have already tested above our mea- 
surements of the vertical velocity dispersions azz 
using a completely different sample of M stars, for 
which the line-of-sight velocities in the direction 
towards the Galactic pole are known. The radial 
velocity dispersions ana can be checked in a sim- 
ilar way. We have determined the first and sec- 
ond moments of the distribution of the pu and pv 
proper motion velocity components of those stars 
in our sample which lie in the cone b > 70°. In the 
upper most counting volume, at 700 pc < Z < 800 



{pjj)2 = 47.9± 0.2 km/s 
= 36.1 ± 0.2 km/s, which 



pc, we find \/ {puPu) 

and V (PvPv) - {pvY 
are in excellent agreement with the result of our 
de-projection formalism (see Table 1). There is no 
independent test of the mixed moment aj^^ within 
our data at this altitude, but we note that the tilt- 
ing angle of the velocity ellips oid of a = 7.3° re- 



ported bv lSiebert et al.l ( 20081 ) implies aj^^ ~ 650 



(km/s)^ , which is even more than what we find. At 
mid altitudes, around 400 pc above the midplane, 
the sample of M stars with known radial velocities 
(described above) provides an independent check 
on the mixed moment crj^^. We have drawn from 
the SpecPhoto database of SDSS DR-7 a set of 
201 stars with b > 70° at heights of 400 < Z < 
600 pc above the midplane, with a mean of (Z) 
= 450 pc. Since the stars lie in the cone b > 70°, 
the line-of-sight velocities refiect their W velocity 
components; we obtain cvi^^ = 29.2 ± 0.7 km/s. 
From the proper motion velocities pu we derive 
ViPuPu) - (pu)^ = 46.8± 1.8 km/s. We have 
then cross-correlated the Vios and pu velocity com- 
ponents, giving {ViosPu) = -275.5 ± 46 (km/s)^, 
which implies a tilting angle of a = 11.2 ± 1.7 
degrees. If the velocity ellipsoid were pointing to- 
wards the Galactic center, one would expect a = 
3.2 degrees. A comparison of these independent 
kinematic data with the seventh and eighth row 
of Table 1 and Fig. [Sj shows excellent agreement 
with the results from the statistical de-projection 
of the proper motion velocities. We conclude from 
this discussion that our measurement of the tilt- 
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ing angle of the velocity dispersion tensor should 
be quite robust. 

Finally, we determine the epicyclic ratio, which 
is related to the slope of the Galactic rotation 
curve as 

dlnvr „ cr^,rf, . ^^^^ 



dlnr 



1 , 



'RR 



Close to the midplane, we find a value of ^j^^^ 
= -0.006 ± 0.016. which indicates a flat Galactic 
rotation curve in the solar neighborhood. Note, 
however, that the quoted error reflects only inter- 
nal errors. Equation (jl7p neglects the asymmetric 
shape of the distribution of the V velocity com- 
ponents of the stars. In order to correct for this 
simplification, third mome nts of the velocity dis- 



tribution must be known ( Cuddeford fc Binnev 



1994I ). which we have not measured here. 



4. Conclusions 

We have analyzed the velocity distribution of 
stars in the Galactic old thin-disk population us- 
ing stellar data from the Seventh Data Release of 
the Sloan Digital Sky Survey. We have constructed 
a sample of about 2 million late- type stars in eight 
cylindrical counting volumes, each of which has a 
height of 100 pc and a diameter of 2 kpc, respec- 
tively, and which were stacked onto one another. 
For each star we have determined a photometric 
distance estimate, and combined this with their 
available proper motions to obtain tangential ve- 
locities. These were then treated as pseudo space 
velocities. We have demonstrated with Monte 
Carlo simulations that our sample provides suf- 
ficient lines of sight through the velocity distribu- 
tion of the stars that the pseudo space velocities 
allow for a precise statistical determination of the 
first and second moments of the three dimensional 
velocity distribution. Close to the Galactic mid- 
plane we find (U) = -8.6 ± 0.2 km/s and (W) = 
-7.1 ± 0.2 km/s, which directly reflect the motion 
of the sun relative to the Local Standard of Rest. 
The mean V velocity component, {V) = -20.0 ± 
0.1 km/s, is much more negative than just the re- 
flex of the solar motion, due to the asymmetric 
shape of the V velocity distribution. Close to the 
midplane, we obtain estimates of the velocity dis- 
persions <7b,r = 32.4 ± 0.1 km/s, cr^^ = 23.0 ± 0.1 
km/s, and azz = 18.1 ± 0.2 km/s, respectively. 
With the exception of (U), all of these results are 



consistent with other determinations, such as anal- 
yses of the Catalogue of Nearb y Stars or Hippar- 
cos stars bv IPehnen fc BinnevI ( 1998f l. However, 
due to the large number of stars we reach an un- 
precedented precision, up to ten times better than 
previously achieved. Moreover, we have a suffi- 
cient sample to trace all quantities as a function 
of height above the Galactic midplane. The ve- 
locity dispersions exhibit a distinct increase with 
height above the midplane. This is a clear indica- 
tion that our sample of M stars contains a mixture 
of young stars with low velocity dispersions and 
old stars with high velocity dispersions, because 
only the latter ascend high above the midplane. 
Since the rotational lag of the stars is tied by the 
asymmetric drift effect to the radial velocity dis- 
persion, -{V) must increase also with height above 
the midplane, precisely as observed. 

We have then studied the orientation of the ve- 
locity ellipsoid in the meridional planes spanned 
by the U, V, and W velocity components. We 
find a vertex deviation of the velocity ellipsoid of 
about 20° to 25°. Close to the Galactic midplane 
the ellipsoid is oriented parallel to the Galactic 
plane. At a height of 750 pc above the midplane 
the velocity ellipsoid is tilted by an angle of about 
a = —20°. This is more than thought previously, 
because a tilting angle of a = —5.4° is expected 
at that altitude, if the velocity ellipsoid points 
straight at the Galactic Center. 

Finally, we determined the epicyclic ratio which 
indicates, close to the midplane, a flat shape of the 
Galactic rotation curve in the solar neighborhood. 

We have refrained in this paper from any de- 
tailed dynamical modelling of our data but hope 
to address this topic in the near future. 
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Appendix A: The Matrices A and B 

The elements of matrix A are given by 

-4(7(7 = 1 — cos^(6) cos^(^) 

Auv = —cos^ {h)cos{l)s\n{l) = Avu 

Auw — — sin(6) cos(6) cos(Z) = Awu 

Aw = 1 - cos^(6)sin^(/) 

Avw = — sin(6) cos(6) sin(/) = Awv 

Aww = cos^(6). (18) 

The elements of matrix B are defined straight- 
forwardly by eq. Illi but are a bit cumbersome 
to write out explicitely. We consider B(nr,nc), 
where nr and nc are the numbers of the rows and 
columns in the matrix B. Next we introduce two 
tables of indices icol and irow and write, if nc 
equal 1, 4, or 6 



B{nr,nc) = —[A{iro'w{nr,l),icol(nc,l)) 



*A{irow{nr,2),icol{nc,2)) (19) 
+A{irow{nr, 3), icol{nc, 3)) 
*A{irow{nr,4:),icol{nc,4))] . 

For the other combinations we have 

B{nr,nc) — —[A{irow{nr,l),'icol{nc,l)) 

'¥A(irow{nr, 2), icol{nc, 2)) 
+A{irow{nr, 3), icol{nc, 3)) 



*A{irow{nr,4:),icol{nc,4:)) (20) 
+A{irow{nr, 1), icol{nc, 5)) 
*A(irow{nr, 2), icol{nc, 6)) 
+A{irow(nr, 3), icolinc, 7)) 
*A{irow{nr,'i),icol{nc,&))] . 



The indices irow and icol of matrix A stand for 
irow, icol — 1, 2, 3, — U, V, W , respectively. 

Appendix B: Error Calculus 

We consider first the moments ([/), {V), and 
(W) . Each proper motion velocity component has 
an individual error en, ey, ew- The error of the 
statistically averaged proper motion velocity com- 
ponents is given by 



N 



2 ^ _ V f 2 

{PU) ,{PV) -.(Pw) J\J2 PU,n,PV,n,PW,n ' 

n=l 



(21) 



where N denotes the number of stars in the count- 
ing volume. We neglect positional errors so that 
matrix elements of (A) have no errors. The errors 
of the first moments of the velocity distribution 
follow then as 



--{U) 



uw 



(pw) 



-i{Avw)-'?eJ 



--{W) 



+ {{A 



WW) 



'{pw) 



The errors of the statistical averages of the combi- 
nations of the proper motion velocity components 
are given by 



N 



-{puPu) 



,2 

{PuPv) 



^ \ " 2 2 

^2 2-^Pu,n^pu,, 



n=l 
N 



1 \ ^ 2 2 I 2 2 

2^'Pv,ri^Pu,^ +Pu,n^pv 



N 

(puPw) ~ ^ ^ PW.n^pu.r, ~^ Pu.n^vw.,1 

n=l 



2 _ ^ „2 2 , „2 2 



-(PvPv) 



E2 2 



(23) 
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Table 3 
Index table irow(nr,j) 



nr 


j=l 


j=2 


j=3 


j=4 


1 


1 


1 


1 


1 


2 


1 


2 


2 


1 


3 


1 


3 


3 


1 


4 


2 


2 


2 


2 


5 


2 


3 


3 


2 


6 


3 


3 


3 


3 



Table 4 
Index table icol(nc,j) 



nc 


j=l 


j=2 


j=3 


,i=4 


j=5 


j=6 


j=7 


j=8 


1 


1 


1 


1 


1 










2 


1 


2 


1 


2 


2 


1 


2 


1 


3 


1 


3 


1 


3 


3 


1 


3 


1 


4 


2 


2 


2 


2 










5 


2 


3 


2 


3 


3 


2 


3 


2 


6 


3 


3 


3 


3 
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n=l 

2 _ ^ \ " 2 2 

n=l 

The errors of the second moments of the velocity 
distribution follow then as 





h2 2 

(puPu) 


f ((S12) 


-1n2 2 

' (puPv) 




U2 2 

' {puPw) 


+ ((S14) 


'"')'4vpv)(24) 


+ m^r 


U2 2 

' {pvPw) 


+ ((Sie) 


/ {pwpwl 




U2 2 

' (puPu) 


f ((S22) 


' (puPv) ■■■ 



and so on. 
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Fig. 2. — Frequency distributions of the proper 
motions in right ascension (top panel) and dech- 
nation (bottom panel) for the stars in each of the 
eight counting volumes. The histograms are each 
normalized to unity and shifted relative to each 
other so that the Z-slice closest to the Galactic 
midplanc appears as lowest and the most distant 
Z-slice at the top. 



Fig. 1. — Difference between the faint and the 
bright normalizatio n of the ab s olute magnitude 
Mr-color relation bv ljuric et al.1 ( 20081 ). as a func- 
tion of (f — i)o color (solid line). Underlaid as white 
histogram is the frequency distribution of the M 
stars in the distance bin closest to the Galactic 
midplanc, while the grey histogram shows the cor- 
responding distribution in the most distant bin. 
Both histograms are normalized to unity (scale on 
the right). 




Fig. 3. — Greyscale diagram of the frequency dis- 
tribution of the reduced proper motions of the se- 
lected M stars versus their (r — z)o colors. The 
distribution is shown on a linear scale, coded from 
black to white. 
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Fig. 6. — Top left panel: The first moments ([/), {V), and (W) as a function of height Z above the 
Galactic midplane (solid dots), as determined in t his work. S olar n eighborhood data (triangles) and the 
results of Dehnen fc BinnevI (1998) (squares) and iVeltz et al.l ( 2008f ) (pentagons) are indicated at Z — 0. 
Bottom left panel: The velocity dispersions c/jfl, a^^, and azz as a function of height Z above the Galactic 
midplane (solid dots) as determined in this work. The solar neighborhood data (triangles) and results of 
IPehnen fc BinnevI ( 1998h (squares) are indicated at Z = 0. Top right panel: Differences of the mean velocities 
([/), (V), and (W) (color coded as in the left panel) determined with the faint no rmalization of the absolute 
magnitudes minus the mean velocities determined with the bright normalization ( Juric et al. 2008h . Bottom 
right panel: Same as upper panel, but for the velocity dispersions crjiji, cr^^, and azz- The right panels 
indicate that the systematic errors of our results are of the order of 2 to 3 km/s. 
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Fig. 4. — Greyscale illustration of the distribution 
of the M stars in our sample as function of height 
above the Galactic midplane and (r — i)o color. 
The density is shown on a linear scale, coded from 
black to white. 
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Fig. 5. — Spatial distribution of the sample stars. 
The upper panel shows the distribution in Galac- 
tocentric radius X and height above the Galactic 
midplane Z. The conical distribution reflects the 
geometry of the survey. The lower panel illustrates 
the projection onto the Galactic plane. The lay- 
ers are shown stacked opaquely onto each other 
viewed from below the Galactic midplane. 
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Fig. 8. — The mixed element —an^ of the velocity 
dispersion tensor as a function of height above the 
Galactic midplane Z. The inset shows the vertex 
deviation of the velocity ellipsoid. 




colors 0.49 < {r — i)o < 0.9. The systematic errors z [pc] 

of our results are of the order of 2 to 3 km/s. 

Fig. 9. — The mixed element — yOfl^z '^^ ^^^'^ veloc- 
ity dispersion tensor as a function of height above 
the Galactic midplane Z. The solid line line is 
a polynomial fit. The inset shows the apparent 
tilting angle of the velocity ellipsoid. 
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